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PREFACE 

This  analysis  was  conducted  for  the  Advanced  Research  Projects 
Agency  as  part  of  RAND's  continuing  study  of  the  propagation  of  laser 
beams  through  the  atmosphere.  It  examines  in  detail  the  characteristics 
of  the  collision-induced  absorption  of  infrared  and  far-infrared 
radiation  by  symnetrical  diatomic  molecules  and  rare-gas  atoms  in 
order  to  provide  a  broader  basis  for  the  theoretical  understanding  of 
collision-induced  absorption,  which  is  an  important  factor  limiting 
transmission  in  some  atmospheric  spectral  windows.  It  is  believed 

that  the  results  reported  here  represent  substantial  progress  toward 
that  end. 
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SUMMARY 

A  quantum  mechanical  calculation  of  collision- induced  absorption 
in  the  infrared  and  far-infrared  spectral  regions  for  mixtures  consisting 
of  symmetrical  diatomic  molecules  and  rare-gas  atoms  is  given.  The  induced 
moment  consists  of  a  short-range  overlap  term  and  a  long-range  quadrupole 
te  Near  molecular  transition  lines  the  quadrupole  term  gives  the 
main  contribution  to  the  absorption,  while  in  the  far  wings  the  overlap 
term  is  dominant.  The  ^"He,  Hj^Ar  and  ^-Ai  systems  are  considered 
and  an  estimate  for  the  absorption  coefficient  near  lQp.  for  N^-N^ 
systems  is  given. 
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I.  INTRODUCTION 

The  purpose  of  this  Memorandum  is  to  calculate  the  binary 
absorption  coefficient  of  infrared  and  far-infrared  radiation  in 
mixtures  consisting  of  symmetrical  diatomic  molecules  and  rare-gas 
atoms.  Atoms  and  molecules  with  no  permanent  dipole  moment,  and 
consequently  no  dipole  absorption  coefficient  proportional  to  the 
density,  are  considered.  The  absorption  mechanism  considered  here 
is  that  due  to  the  dipole  moment  induced  between  a  pair  of  colliding 
atoms  or  molecules.  For  small- impact  parameters  the  induced  moment 
arises  as  a  consequence  of  the  distortion  of  the  electron  cloud  during 
the  collision  of  two  particles,  atoms,  or  molecules,  while  for  larger- 
impact  parameters  the  induced  moment  arises  from  the  polarization  of 
one  particle  by  the  quadruple  field  of  the  other.  The  induced  dipole 
moment  is  expected  to  be  nonzero  when  the  colliding  pair  are  within 
a  few  molecular  radii  of  each  other.  Thus,  absorption  takes  place 
within  the  duration  of  a  collision  and  the  probability  of  absorption 
is  proportional  to  the  product  of  the  densities  of  the  two  species. 

Collision-induced  absorption  in  rare-gas  mixtures  was  discovered 
experimentally  by  Kiss  and  Welsh^  and  further  studied  by  Bosomworth 

/n\ 

and  Gush.  Recently,  Levine  and  Birnbaum( 3,4)  have  made  a  theoretical 
calculation  of  induced  absorption  in  rare-gas  mixtures  using  a  simple 
model  to  describe  the  variation  of  the  collision-induced  dipole  moment 
with  internuclear  separation.  They  obtain  an  analytic  expression  for 
the  absorption  spectrum  which  agrees  well  with  the  experiments. 

Previously,  J.  Van  Kranendonk^  ^  has  presented  a  theory  for  the  integrated 
absorption  coefficient  in  diatomic  molecules.  Interest  here  is  in 
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the  ‘.bsorption  coefficient  as  a  function  of  the  frequency  of  the  incident 
radiation  field  and  not  in  this  function  integrated  over  all  frequencies. 

The  energy  absorbed  during  a  col  .stion  between  two  dissimilar 
rare-gas  atoms  goes  into  increasing  the  translational  kinetic  energy 
only.  The  effects  on  the  induced  absorption  due  to  the  rotational 
transitions  that  occur  in  the  diatomic  molecule  are  investigated  here. 

The  frequency  of  the  applied  radiation  field  is  taken  to  be  less  than  the 
frequency  of  the  first  excited  vibrational  level  of  the  diatomic 
molecule  (4159.4  cm  *  for  and  2330.7  cm  *  for  see  Ref.  6); 
consequently,  vibrational  effects  on  the  absorption  need  not  be  con¬ 
sidered.  The  absorbed  energy  in  this  c'se  goes  into  both  rotational 
and  translational  energy. 

In  Section  II  a  formal  expression  for  the  absorption  coefficient 
is  written  and  the  phenomenological  induced  dipole  moment  is  introduced 
and  discussed.  In  Section  III  an  analytic  expression  for  the  absorption 
coefficient  is  derived,  and  in  Section  IV  an  estimate  is  made  of  the 
parameters  for  H^-rare  atom,  N^-rare  atom,  and  collisions,  and 

numerical  values  of  the  absorption  coefficient  are  given  as  calculated 
as  a  function  of  frequency  for  these  cases. 

These  calculations  are  aimed  at  providing  a  broader  basis  for  the 
theoretical  understanding  of  collision-induced  absorption,  which  is  an 
important  factor  limiting  transmission  in  some  atmospheric  spectral 
windows.  For  example,  collision-induced  absorption  in  0 ^  is  dominant 
at  the  Nd  laser  wavelength  of  1.06u,  and  collision-induced  absorption 
in  the  far  wings  of  HgO  lines  is  important  at  the  CO 2  laser  wavelength 
of  10. 6y..  These  absorptions  involve  more  complicated  mechanisms  and 
require  further  development  for  their  treatment. 
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II.  FORMAL  EXPRSSSIC,  FOR  THE  INDUCED  ABSORPTION  COEFFICIENT 

2 

Gaussian  units  are  used  with  ft  ■  c  «  1;  then  e  2!  1/137.  The 
absorption  coefficient  cr  =  (2TT/c)r,  where  F  is  the  transition  ratt  . 

That  is, 


a  (as) 


'  *  *>if|2<pi‘pf>  6(<Dfl-co^ 


(1) 


In  Eq.  (1),  V  is  the  volume  of  the  sample,  £  is  the  induced  dipole 

moment  of  the  system  of  particles,  e  is  the  unit  polarization  vector  of 

the  incident  radiation  field  of  frequency  cot  *  ®f”Pi  *-8  ener8y 

difference  between  a  final  state  f  and  an  initial  state  i  of  the  molecular 

system  respectively.  P^  and  are  the  statistical  weights  of  the  initial 

and  final  particle  states.  Also,  Eq.  (1)  indicates  a  sum  over  final  states 

(denoted  by  £)  and  an  average  over  initial  states  (denoted  by  . 

xne  appearance  or  the  factor  (P  -Pp  takes  into  account  both  absorption 

(the  term)  and  the  reverse  process,  spontaneous  emission  (the  P^  term). 

It  is  assumed  throughout  that  the  temperature  T  of  the  system  is  so  low 

that  all  of  the  particles  may  be  taken  to  be  in  their  electronic  ground  states. 

Consider  a  system  of  rare-gas  atoms  (density  n^)  and  symmetrical 

diatomic  molecules  (density  nj).  The  dipole  moment  arises  from  binary 

*★ 

collisions  between  atoms  and  molecules  and  between  molecules.  Con¬ 
sider  first  the  term  in  or  proportional  to  nj^,  that  is,  absorption 


Explicit  expressions  for  P.  and  F,  are  given  in  Eq.  (7)  below. 

^  ^  * 

There  can  be  no  net  induced  dipole  moment  for  two  identical  rare- 

gas  atoms  (in  the  same  electronic  state)  since  the  system  is  completely 
symmetrical  about  the  center  of  mass.  However,  this  is  not  the  case  for 
two  diatomic  molecules. 


■Jiiuj|i;l|iudjlini<illiiii!iiiiiii 
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by  l.hii  induced  moments  arising  from  the  collision  of  a  rare-gas  atom 

and  a  symmetrical  diatomic  molecule.  Later  on,  the  term  proportional 
2 

to  n  (i.e.,  absorption  arising  from  binary  molecular  collisions)  i" 


estimated.  It  is  assumed  that 


Vnx  Vn2 


^V-’V  =  Z  Z  ^(Rij^ 

i=l  j=2 


i.e.,  the  dipole  moment  arises  from  binary  collisions  and  is  a  function 

of  the  distance  between  the  centers  of  mass  (R^)  the  two  colliding 

particles.  The  wave  functions  describing  the  centers  c"  mass  motion 

of  the  particles  can,  in  the  first  approximation,  be  taken  as  plane 

waves.  Then  by  a  straightforward  procedure  Eq.  (1)  may  be  reduced  to 

(4) 

a  form  expressing  a  in  terms  of  two-particle  matrix  elements.  '  Thus 


cr  =  4tt 


»l'2»(  Y  !(?(*) 'e)ifJ2  (Pi-Pf)  6(u>fi-«))\  (3) 

\f  /  1 


where  now  the  initial  and  final  states  refer  to  the  two-particle  system 
interacting  with  the  radiation  field, 


£(R)*e|  u  ■  I  cpf(R)V(R)-e)cpi(R) 


where 


^(R)  -  eik-R  u^R)  +  (0)  =  eik‘R  *^(0) 


cp  yrO  -  eik/*Ruk/(R)  ^(fl)  2:  ejk,#R  *^(0) 

Experimentally  these  two  cases  may  be  distinguished  by  .or  example) 
arranging  n  »  so  that  the  absorption  arises  predominately  from 
collisions  Between  diatomic  molecules  and  rare-gas  atoms. 

The  *i(R)  terra  in  the  relative  two-particle  wave  function,  which  is 
determined  by  solving  the  Schoedinger  equation  with  the  two-particle  inter¬ 
action,  is  taken  to  be  unity  in  the  first  approximation.  This  corresponds 
to  taking  straight  line  collision  paths  in  the  classical  approach. v3>4) 
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tj  are  the  rotational  wave  functions  and  k  and  k*  are  the  relative 
momenta  of  the  colliding  particles  before  and  after  the  collision, 
respectively. 

For  the  atom-tnolecule  system, 


o 

E.  -  k  /2m  +  E 

Ji 

E ,  -  k,2/2m  +  F, 

► 

m  ■  reduced  mass  «  m]m2/(mi+ra2^ 


Ej  *  rotational  energy  ■  BJ(J+1) 


(6) 


where  J  is  the  rotational  quantum  number  and  B  is  the  rotational  constant 

of  the  molecule.  The  statistical  weights  of  the  initial  and  final  particle  states 

are  given  by 


and 


-E  /hT  . 

Pt  -  (2J1+l)gi  e  1  (l/2nnmT)J// 

"Ef/HT  1/2 

Pf  -  (2Jf-*-l)gf  e  (l/2nmHT)J/2 


(7) 


where  k  is  the  Boltzmann  constant,  and  g  and  g^  are  the  statistical  weights 
for  the  nuclear  spin  states  of  the  diatomic  molecule.  It  is  well  known 
that  the  relative  statistical  weights  for  the  symmetric  and  anti¬ 
symmetric  spin  states  are^ 


8(sym) 


/g 


(antisym) 


(I+l)/l  for  I  integral 
i/l+l  for  I  half-integral 


(8) 


where  I  is  the  nuclear  spin.  Upon  combining  Eqs.  (3) -(8), 


m 
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«(<»)  -  AtTZn1n2a)(l/2nmHT)3/,2(2J1+l)"1  V  P  jiV.  P  d3k/ 

(2TT>J  <2TT)3 


*  I 


,-(lc  /2mnT) 


- —  f(2J  +l)g 

J,(J.+1)  <•  1  1 


f  1  r  ,  ■  S  JiWi+D 

£  (2^+1) 8l  e  mT  1  1 

i 

-  <yf+l)gfe-“*T}  *  J‘U‘+1)  h„ 


(9) 


where 


“  '  ®fi  *  «>  '  S  (k'2_k2)  '  E 


JfJi 


(10) 


and 


In  Eq.  (9)  the  summations  over  the  initial  and  final  continuum  states 
k  and  k7  have  been  replaced  by  integrations  over  d3k  and  d3k' 
respectively. 

To  proceed,  an  analytic  form  is  specified  for  the  induced  dipole 
moment.  The  induced  dipole  moment  of  a  pair  of  colliding  particles  may 
be  resolved  into  two  parts.  One  part  is  the  overlap  dipole  moment 
which  results  from  the  distortion  of  the  electron  cloud  during  close 
collisions  (r  <  collision  diameter).  The  other  part  of  the  induced 
dipole  moment  arises  from  the  polarization  of  one  particle  by  the 
quadrupole  field  of  the  other  and  vice  versa.  The  quedr  jole-induced 
dipole  moment  arises  from  long-range  collision  (r  i  collision  diameter). 
These  two  parts  of  the  moment  are  discussed  separately. 
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(3  4) 

In  view  of  the  close  agreement  of  Levine  and  Birnbaum's  * 

theory  with  the  experimental  results,  we  use  the  same  functional 

form  for  the  overlap- induced  dipole  moment  between  the  rare-gas  atom 

and  each  nucleus  of  the  diatomic  molecule.  The  form  they  use  is 
2  2  -1 

ji(r)  «=  ^cYr  exP  (“Y  r  )•  w^ere  Y  an<*  are  t*ie  range  and  strength 
respectively  of  the  dipole  moment.  In  the  cares  considered  here 

°_  i 

Y  ~  (1-2)  A  .  That  is,  the  important  values  of  r  are  of  the  order 
of  the  internuclear  separation  of  the  me  ecule.  For  distance  in 
this  range  the  rare-gas  atom  "sees"  each  nucleus  separately.  This, 
however,  would  not  be  the  case  when  y  ^  is  much  larger  than  the  inter- 
nuclear  separation. 

In  the  present  calculation  any  interparticle  potential  is  ignored, 
and  hence  straight  line  collision  paths  are  assumed.  Realistically, 
a  potential  exists  between  the  particles  which  contains  a  repulsive 
part  that  prevents  the  two  particles  from  approaching  too  close  to 
one  another,  so  there  is  actually  a  small  moment  for  distances  within 
the  range  of  the  repulsive  part  of  the  interparticle  potential.  The 
factor  (yr)  in  the  moment  cuts  off  che  moment  for  r  $  y  *.  Similarly, 
for  large  interparticle  separation,  when  the  electron  clouds  do  not  over¬ 
lap  there  is  very  little  contribution  to  the  overlap  moment.  The  factor 

2  2  >  -1 
exp  (-y  r  )  cuts  off  the  moment  for  r  ~  y  .  Strictly  speaking,  the 

repulsive  and  overlap  ranges  are  not  equal.  In  this  calculation,  as 

in  Refs.  3  and  4,  a  single  range  is  used  to  specify  the  induced  moment. 

By  using  the  modified  gaussian  model  for  the  induced  overlap 
moment  Levine  and  Bimbaum  were  able  to  obtain  a  simple  analytic 
expression  for  the  absorption  spectrum.  By  matching  the  theoretical 
and  experimental  results  (for  Ke-Ne,  He-Ar  and  Ne-ar  systems)  at  the 
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peak  of  the  spectrum  they  obtained  values  of  y  and  Using  these 

values  of  v  and  they  obtained  good  agreement  over  nearly  the  entire 
1  o 

range  of  the  spectrum  (especially  for  He-Ar  systems).  For  the  cases 
considered  here,  H  and  N  /rare-gas  atom  absorption,  no  experimental 
results  exist  as  yet.  However,  using  the  mathematical  form  for  the 
moment  that  was  used  by  Levine  and  Birnbaum  leads  to  analytic  expres¬ 
sions  for  the  absorption  spectrum  in  terms  of  and  y.  Then  physical 
arguments  are  used  to  estimate  values  for  and  y  in  terms  of  the 
values  obtained  by  Levine  and  Birnbaum  for  the  rare-gas  cases.  In 
view  of  their  good  agreement  with  the  experimental  results  this 
procedure  seems  justifiable. 

The  main  contribution  to  the  induced  dipole  moment  arises  for 

o-l 

values  of  r  where  ry  ~  1.  In  the  cases  considered  here,  y  ~  (1-2)A  . 

That  is,  the  important  values  of  r  are  of  the  order  of  the  size  of 
the  atom  or  molecule  (the  electron  clouds  overlap  and  hence  distort 
the  charge  distri>ution) . 

The  quadrupole-induced  dipole  moment  is  equal  to  the  quadrupole 
field  of  particle  one  at  the  position  of  particle  two  times  the 
polarizability  of  particle  two,  and  vice  versa.  In  the  case  of  a 
diatomic  molecule/rare-gas  atom  collision  there  is  only  one  terra, 
since  rare-gas  atoms  possess  no  quadrupole  moment.  Qualitatively 
the  overlap- induced  moment  contributes  t>.  the  high-frequency  part  of 
the  absorption  spectrum,  since  the  overlap  moment  arises  from  close 
collision.  On  the  other  hand  the  quadrupole-induced  moment  contributes 
to  low  frequencies  in  the  Q  and  0  branches  and  frequencies  near  the 
rotational  transitions  in  the  S  branches,  since  the  quadrupole-induced 
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moment  is  due  to  distant  collisions.  (Near  the  rotational  transition 
frequencies  the  momentum  transfer  is  small.) 
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III.  CALCULATION  OF  THE  ABSORPTION  SPECTRUM 

OVERLAP- INDUCED  DIPOLE  MOMENT 

The  induced  dipole  moment  between  a  colliding  rare-gas  atom  and 
a  symmetrical  diatomic  molecule  is  taken  as  a  superposition  of  the 
atom  interacting  with  each  of  the  two  atomic  components  of  the  molecule 
respectively.  That  is, 

£<o)(r)  -  ^(rp  +  U(0)(r2)  (12) 

where  r^  and  r^  are  the  separations  of  the  atom  from  each  of  the  two 
nuclei  respectively.  Choosing  a  coordinate  system  with  the  z-axis 
parallel  to  the  polarization  vector  of  the  incident  radiation  field  and 
calling  the  molecular  internuclear  separation  vector  2a  gives 

.  2r2  .  2r2 

u^(r)*e  ■  uQv j^ri>  e  1  +  e  2]  (13) 

"  z  z 

where 

-»  _•  -*  .  -*  .  -» 

r^  -  r  -  a  and  r^  *  r  +  a 

First  the  matrix  element  of  u«e  between  the  initial  and  final 
continuum  states  is  calculated.  Suppressing  the  rotational  quantum  numbers 
for  brevity,  we  have 


(?(0)-eVk  =  Jd3r  e'ii?''Ve)  e1*7 
=  J  d3r  (?-e)  e'^‘r 


(14) 


where 


q  -  1'  -  TZ 
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Substituting  Eq.  (13)  into  Eq.  (14)  gives 

<u  -S)k'k  "  <uzVk  "  |  d3rj.  ^  1  «,  e  V  1 


+  e"1(J*a  T  H3  1?'?2 
+  C  |4r2*  *2*  2} 


-  vu0  .*S  J  d3r  .  ^  e-V2r2  +  teraB  ;  „  -  (u) 

The  integral  appearing  in  Eq.  (15)  is  elementary  and  easily  found  to  be:* 

f  d3r  .  e^'?  e-v2*2  .  Hi  -5  -«W) 

J  2i  T  xz 

Substituting  this  value  into  Eq.  (15)  gives 


<*°-s>k'k  ■  d-(Qi/4/> 


2/a.  2. 


(16) 


+  terms  a  -»  -a 


Next  one  must  calculate  the  matrix  element  of  Eq.  (16)  between 

t/i  and  final  rotational  states  of  the  molecule.  Suppressing  the 
subscripts  and  k, 


<»!> 


JA 


.,A 


where  the  initial  and  final  rotational  states  are  referred  to  the  z-axis. 
Suppressing,  for  brevity,  terms  that  do  not  depend  on  angular  coordinates, 

* 

Due  to  the  rapid  convergence  of  the  integrand  for  r  »  v"1  nealieible 
error  is  introduced  upon  extending  the  upper  limit  of  the  r  integration 


I 
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we  find 


\ 


g^*®! terras  a  -•  -a  (17) 


For  frequencies  of  the  order  of  the  peak  in  the  absorption  spectrum, 
i  i  <  iQ.a 

|Qa|  ~  1  and  hence  one  may  expand  e  and  retain  the  first  nonvanishing 

term.  In  the  case  of  diatomic  molecule/ rare-gas  atom  collisions  the 

i(J*  a 

main  contribution  comes  fr-r:..  setting  e  al,  That  is,  for  diatomic 
molecule /rare-gas  atom  collisions  the  overlap- induced  moment  is 
independent  of  the  orientation  of  the  molecule.  On  the  other  hand, 
for  molecular-molecular  collisions  involving  identical  molecules  (in 
the  ground  electronic  and  vibrational  states)  the  matrix  elements  of 
the  angle- independent  part  of  the  overlap- induced  moment  are  zero,  since 
the  colliding  system  possesses  a  center  of  symmetry.  Hence,  the  lowest- 
order  nonvanishing  term  for  systems  such  as  Nj-Nj,  Hg-Hg  comea  ^rom 
the  (Q'a)  term  in  the  expansion  of  e^  &  (see  Eq.  (17)).  The  angle- 
dependent  part  of  the  overlap  moment  is  small, ^  and  hence  for  fre¬ 
quencies  of  the  order  of  the  peak  the  overlap  moment  contributes  to 

the  Q  branch  only.  For  frequencies  that  |Qa|  ~  1  (tail  of  the 

,  ■ — * 

spectrum),  a  power  series  expansion  of  e  *a  is  not  valid  and  the  over¬ 
lap  moment  is  strongly  angle-dependent.  For  example,  in  N  -N  absorP“ 
tion  the  peak  in  the  spectrum  is  approximately  100  cm  *  and  arises 
from  quadrupole-induced  dipole  moments.  Near  1000  cm  ^  (~  (#2  laser 
wave  number)  where  Qa  >  1,  the  main  contribution  to  the  absorption  is 
expected  to  come  from  the  angle-dependent  terms  in  the  overlap-incb'ced 
moment.  For  the  present,  attention  is  focused  on  diatomic  molecule/rare- 
gas  atom  collisions  near  the  peak  of  the  spectrum.  Later  on,  the 
absorption  in  N„-N0  collisions  near  the  C0o  laser  wavelength  is  estimated. 
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— *  — * 

Setting  e  =“  1  in  Eq.  (16)  ,  the  overlap  matrix  element  is 


given 


by 


2  /  .  2N 


/4v  >8j£J15Mf«1  (18) 


QUADRUPOLE- INDUCED  DIPOLE  MOMENT 

For  diatomic  molecule/rare-gas  atom  systems  the  long  range 
quadrupole- induced  moment  ,  which  can  be  inferred  from  Ref,  7, 
is  given  by 


t.-(q)  „  lVo  4TT. 

r4  {  5  ' 


m*+2 


I  W  Wnr> 


m*-2 


s*er  for  r  £  o  (19) 


where  the  z-axis  is  taken  parallel  to  e,  e  is  a  unit  vector  in  the 

r 

r-direction,  and  0^  are  the  angular  coordinates  of  r  and  a 

respectively,  is  the  polarizability  of  the  rare-gas  atom,  and  q 
is  the  quadrupole  moment  of  the  diatomic  molecule.  The  assumption  is  made 
that  the  quadrupole-induced  moment  has  meaning  only  for  values  of  r  £  a, 
which  is  the  collision  diameter  of  the  colliding  particles. 

Upon  taking  matrix  elements  of  Eq.  (19)  between  initial  and  final 
continuum  states, 


where 


12tt 


a  q 
oMo 


Y20<na> 


e*I 


(20) 


I 


*20<n>  '1Q'r  V 

r 


(21) 


and  Q  ■  k  -k.  I  is  a  vector  and  must  be  proportional  to  Q,  That  is, 
I»IQ  where  I-((J«f)/Q2.  Then 


3*?  -  J  y20(C1) 
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Taking  the  z-axis  of  integration  along  Q  and  expanding  e^’r  in  terms 
of  spherical  harmonics, 


-  -  dr  BVQr) 

q-i-  iv«;  q  f 


(22) 


where  *-s  t^e  spherical  Bessel  function  of  second  order.  The  lower 
limit  on  the  r  integration  has  been  set  equal  to  the  collision  diameter 


a  appearing  in  the  interparticle  potential. 

Strictly  speaking,  setting  the  lower  limit  arbitrarily  equal  to  a 

is  incorrect.  Upon  introducing  the  two-particle  distribution  function 

into  the  Integrand  of  Eq.  (22),  one  can  now  integrate  over  all  r.  The 

two-particle  distribution  function,  in  the  low  density  limit  and  at 

-Vfr) /rT 

temperatures  considered  here,  is  proportional  to  e  v  '  where  V(r) 
is  the  interparticle  potential.  Any  realistic  interparticle  potential 
(for  example  the  Lennard-Jouss  potential)  becomes  very  large  for  small 
r  (r  «  size  of  molecule),  and  hence  (e  '  '  «  1) ,  a  negligible 

contribution  to  the  integral,  come3  from  small  r.  For  r  ~  size  of  the 
molecule,  V(r)  **.  0  (i.e,,  the  two-particle  distribution  function  is 
constant).  The  main  contribution  to  the  integral  comes  from  r  “«a. 

In  view  of  other  approximations  in%'olved  in  the  calculation  it  is  felt 
that  unnecessary  mathematical  complexity  is  introduced  by  using  the 
exact  two-particle  distribution  function.  Therefore  we  use,  as  an 
approximation,  a  distribution  which  is  zero  for  r  <  o  and  is  equal  to 
unity  for  r  2  o. 

For  frequencies  tu  ~  *T  for  the  Q  branch,  and  frequencies  near 
the  rotational  transition  frequencies  for  the  S  branch,  Q,  the  momentum 
transfer,  ~^/2m>tT.  For  rT  approximately  equal  to  room  temperature  and 
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v&lues  of  d  and  c  considered  here,  Qa  »  1  and  hence  the  asymptotic 
value  of  J2  may  be  used.  Setting  x  -  Qr,  Eq.  (22)  becomes 


-*  —  2  dCj2(x)] 

Q»  I  “  v20tt  iQ  f  - - 

■I  rv 


Qg  x 


2:  -^/20tt  i  Q2  ^-n  +  0(Qc)‘4  +  ... 


Therefore, 


I*e  -  (Q*e)(Q*I)/Q^ 


m 


«otoi^3/20i-t)  !iS-f  Y  (0)  (24) 

v  (Qct) 

Taking  matrix  elements  with  respect  to  the  rotations]  wave  functions 
(suppressing  factor,  that  do  not  depend  on  rotational  quantum  numbers), 


it  is  found  that 


(Vf  l^<q)  K*<i)~  (•>£*£  |K2oIJA) 


\5(2J^+1)  "H 

4tt(2J,+1)  c(J12Jf>Mi0)c(Ji2jf, 
-  *■  J 


00)  (25) 


where  the  c's  are  the  well-known  clebsch-Gordan  coefficients. (8)  Fran 
Eq.  (25)  the  selection  rules  for  quadrupole  absorption  arc  easily 
obtained]  they  are  4J  .  -  0,  *  2  and  ^  =  0.  Fro.  E„.  (24) 


and  (25)  it  is  found  that 
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12nia  q  (CJ*e) 
o  o 


sin  Qg  (2Ji+1\* 

(Qc)3[2Jf+1J 


X  C(J12Jf,Mi0)C(Ji2Jf,00)  (26) 


where  Q  «  k'  -  k. 

CALCULATION  OF  THE  ABSORPTION  SPECTRUM 

According  to  Eq.  (5),  It  Is  necessary  to  square  the  matrix  elements  of 
U  =  jj.e  (n  ■  and  then  sum  over  final  states  (rotational  and 

continuum)  and  average  over  Initial  states  (rotational  and  continuum). 
Consider  the  rotational  part  first:  suppressing  notation  that  does 
not  depend  on  rotational  quantum  numbers,  the  average  for  a  given 
initial  rotational  state  takes  the  form 

U2M2Ji+l)'i  y  tt-e  IJ^y)!2 

*VMf 

where  n*e  ■  (g/  u^)*e.  From  Eqs.  (18)  and  (26)  it  is  found  that 


U2  -  (Q-l)2!*2  e-<oW),  +  A2(i>£)  \ 

l  0  JfJi  q  (Qa)CJ 


—  (Q-e)2  (R0+Rq) 


(27) 


where 


-3/2  -4 
u  n  y 


and 


A2(i,«  -  (12™o,o)2  C2(Jl2J.,00) 


No  cross  term  between  overlap  and  quadruple  moment  appears 
in  Eq.  (27).  The  cross  term  is  proportional  to 

J  C(J12Jf,.V)5,  ,  -  7  W^.M.C)  -  0 


Ml-Mf 


'iWf 


M. 
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The  first  te  m  in  Eq.  (27)  is  the  contribution  from  the  overlap  term, 
while  the  second  is  the  contribution  from  the  quadrupole  term.  The 
angle- independent  overlap  term  contributes  to  the  Q  branches  only, 
while  the  quadrupole  term  contributes  to  the  S,  Q  and  0  branches. 

To  evaluate  the  integrals  over  £  and  k'  indicated  in  Eq.  (9), 
variables  are  changed  from  k,  ic'  to  k,  Q  respectively.  Then  (suppres¬ 
sing  factors  that  do  not  depend  on  k  or  Q) , 

*(«#)  ~  j*  -^3  J  -^3  e"€k  U2  6f(2m)'1(Q2+2Q-k-A2)l  (2 

J(2n)3J(2n)3  L  J 


where 


A  =  2m(u<  -  E  ) 
Jfi 


l/2mHT 


Holding  Q  fixed  for  the  k  integration, 


2  2  2 

Of(m)  ~  (2rr)  2(m/Q)  J  k  dkd(cos  y)e",k  6(cos  x+  2qX') 

where  x  is  the  angle  between  k  and  Q  and  the  Q  integral  is  suppressed. 
Performing  the  integral  over  cos  x, 

00  2 

<*(  )  ~  (2-)_2(m/Q)  f  k  dk  e"*k  (2 


where 

k  = 

kQ  2Q 

The  limits  on  the  k  integral  are  determined  by  requiring  |cos  y|  £  1 

and  k  1  0.  The  integral  appearing  in  Eq.  (29)  is  elementary,  and 

-2  v2 

or(«))  ~  (2tt)  (m/2Qt)e"*  Q. 
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Hence  from  Eq.  (28), 


,  .  .3-  -ek  2 

a(u>)  -  (2n)0(m/2c)  J  ^  e  *  (Q-?)  (Ro+Rq) 

R  and  R  are  defined  in  Eq.  (27).  When  the  z-axis  of  integration  is 
o  q 

chosen  to  be  parallel  to  e,  the  angular  integration  becomes  elementary 

* 

and  gives 


®(<o) 


CO  ■» 

-  (2TT)"4(m/3e)  f  Q3  dQe  Q  (R  +R  ) 
!0  °  1 


2  .  -CQ2C— y  +  f)  + 

(2H)-4  (W3.)  /2  f  Q3  dQ  {A2  .  ^ 

J0  *■  ° 


2  -4+4 

+  A2  SiS-Sa  .  4Q2 

’  (Qa)6 


} 


(30) 


12tt)-4  (mfr/3e)  eeA  /2  f  x3dx  /a2 

'o  »■  ° 


-[x2(a+b)  +  -^3 
x 


A2  sin2Bx  'b(x  +  2^ 

+  -9 - e 

(rx)6 


} 


2  2  2 

where  a  »  A  /2y  ,  b  “  eA  /4,  and  C  ■  Act.  Two  types  of  integrals  appear 
in  Eq.  (30).  The  first  is  of  the  form 


Strictly  speaking,  the  upper  limit  on  the  quadrupole  term  should  be 
<*raax  ~  Ztt'0*  However,  due  to  the  rapid  convergence  of  the  integrund  for 

nnor  1  _ a. 


large  Q  the  upper  limit  has  been  set  to  ». 
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I,  - 


r  x3  dx 

'o 


-(ft  X2  + 


“  -(«y  +  ~) 


i  f  y  dy  e 
J0 


.^{r-^y+?)} 

M  l',0  y  J 

a2Ko(2  Jrf)  * 

— — 

-  (p/a)  1^(2  yjsj 


where  Kn  Is  the  modified  Bessel  function  of  the  second  kind  of  order  n. 

The  second  integral  appearing  in  Eq.  (30)  cannot  be  expressed  in 
closed  form.  However,  the  limiting  cases  of  large  and  small  C  lead 
to  closed  form  expressions. 

„  “(ox2  +  \) 

h  ”  f  e  X  sin2  01 

•  0  x 

2 

For  C  «  1  the  sin  term  varies  slowly  over  the  region  where  the 
exponential  term  is  nonzero.  Applying  the  method  of  steepest  descent 
to  I2  it  is  found  that 

I2  “  (ft/3)^  sin2[c(S/ft)^]  Kq(2  -/ffP),  C  «  1 


a  C2Kq(2  v^Pj 


2 

For  C  »  1  the  sin  term  oscillates  rapidly  over  the  *-egion  where  the 
exponential  term  is  nonzero,  and  hence  its  average  value,  may  be 
substituted,  yielding 


See  Ref.  9. 
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I2  k  1^(2  v^B),  C  »  1 


The  integrals  appearing  in  Eq.  (30)  may  be  obtained  by  substituting 
a  -  afb,  g  ■  b  in  1^  and  or  ■  B  ■  b  in  I The  result  is 

2  A? 

«(o>)  -  (2rr)"4(tnA4/3e)  /2  {~  ^(21^1+6)  +  ~f  I.J  (31) 

C 

where 


6  -  a/b  *  4mKT/y 


and 


.  -(*2+*^)b 

*^X  e  X  sin2  Cx 


j*  dx 

Jn  3 

”0  X 


-*  C2KQ(2b)  for  C  «  1 
-  h  Kx(2b)  for  C  »  1 


Replacing  all  of  the  factors  previously  withheld,  it  is  found 


that 


3 

n^n^m  <u 


«(cd)  -  r  7  {(““Ej  j  ) 

3ttJ(2h»iwT)*Z  .  {  f  i 


1  e  2kT 


3tt(2h»iwT)^Z  T  % 
Ji,Jf 


(2Jf+l)  eWHrj  e  KT 


^[(2^+1) 


Ji(Ji+1)J 


y  (32) 


p"  K  (2b^+6)  e 
Ll+6  1  Jii 


,6  q 


]} 
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where 

b  -  (»-E.  T  )/4kT 
2 

6  ■  4mHT/*y 
C  ■  bo 

Z  -  £  <2Ji+l)gi  e"  RT 
Ji 

Ao’  Aq  and  \  are  «iven  in  Eqs.  (27)  and  (31)  respectively.  In 
obtaining  Eq.  (32)  has  been  set  equal  to  g^,  since  the  interaction 
does  not  change  the  nuclear  spin. 

Examining  the  overlap  term  in  Eq.  (32)  it  is  found  that  it  is 

identical  to  the  results  given  by  Levine. (4)  This  is  expected,  since 

in  the  present  approximation  the  molecular  overlap  term  is  independent 

of  molecular  orientations  and  hence  is  formally  the  same  as  the  rare- 

gas  case.  Examining  this  term  further  it  is  found  that  the  peak  u) 

o 

in  spectrum  is  given  by  ^  -  6hT[<1+6)*-13_1  where  6  „  Um^T/y2.  For 
low  temperatures  (kT  «  *To  -  Y2/4m)  cr  is  independent  of  temperature, 
while  for  high  temperatures  (*T  »  hTq  =  Y2/4m)  a  -  (xT)*  exp  [-(u^/yCnT)*)  ]. 
The  values  of  xTq  for  the  systems  considered  here  are  100  cm"1,  40  cm”1, 
and  -  6  cm  for  H^He,  H^Ar,  and  N^Ar  respectively.  Room  tempera¬ 
ture  (xT  =»  200  cm  )  lies  in  the  high  temperature  region  for  the 
cases  considered  here.  Also,  the  temperature  dependence  of  the 
location  of  the  peak  in  the  absorption  spectrum  is  given  by 
(kT  «  hTq)  0)°  ~  constant  independent  of  temperature  for  low  tempera¬ 
tures,  and  (xT  »  hTq)ujo  —  (kT)^  for  high  temperatures.  That  is, 
for  low  temperatures  the  peak  in  the  spectrum  remains  fixed  as  the 
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temperature  is  varied,  while  for  high  temperature  the  peak  moves 
toward  higher  wave  numbers  as  the  temperature  is  increased  (and  vice 
ve  rsa) . 

The  width  of  the  absorption  spectrum  Atu  is  approximately 
k  - 1 

6hT[(1  +  1/Tq)*  -  1]  .  For  low  temperatures  (T  «  Tq)  ,  Au)  ~  12*To 

(i.e.,  independent  of  temperature),  while  for  high  temperatures 

l  k 

(T  »  T  )  ,  Atu  ~  6(kTq)^  (xT)^  (i.e.,  proportional  to  the  square  root 

of  the  temperature).  These  two  limiting  cases  may  also  be  obtained 

directly  from  physical  arguments.  The  low  temperature  case,  T  «  T^, 

implies  P  «  y  where  P  is  the  relative  momentum  of  the  colliding  particles. 

The  induced  moment  gives  rise  to  an  attractive  force  which  increases  the 

2 

relative  kinetic  energy  from  a  very  small  value  (P  /2m)  at  large 

2  -1 

distances  to  a  value  AE  ~  y  /2m  at  distances  of  the  or C-."  y  .  From 

the  uncertainty  principle  (AEAt  ~  1),  the  transition  rate  (^  At  *) 

-1  2 

and  hence  the  width  of  the  spectrum  is  Acu  =  2nAE  =  At  ~  TTy  /m  =  4tthT 

1  o 

(i.e.,  independent  of  temperature).  For  the  high  temperature  case, 

T  »  Tq,  we  have  Au:  ~  2r?At  In  this  case  the  rate  is  approximately 

\  "1 

(nT/m)  and  the  distance  approximately  y  (as  can  be  inferred  from 

J.  it 

the  uncertainty  principle).  Therefore,  Ad)  ~  4n(xT)  *  (hTq)  . 

Numerical  estimates  of  the  overlap-induced  dipole  parameters  in 
^2’  ^/rare-gas  atom  systems  and  quantitative  results  for  the  absorp¬ 
tion  coefficients  for  these  cases  are  given  in  the  next  section. 
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IV.  ABSORPTION  SPECTRUM  OF  HARARE  ATOM  AND  N^-Ar  MIXTURES 

OVERLAP  DIPOLE  PARAMETERS  FOR  DIATOMIC  MOLECULE/RARE-GAS  MIXTURES 

In  order  to  obtain  numerical  results,  estimates  of  the  phenome- 

(3  4) 

nological  parameters  p  and  y  are  required.  Le /ine  and  Bimbaum  *  ' 
have  determined  these  parameters  for  three  cases  of  rare-gas  atoms; 
values  are  given  in  the  following  table. 


DIPOLE  PARAMETERS  FOR  RARE-GAS  MIXTURES 


Parameter 

System 

He-Ne 

He-Ar 

Ne-Ar 

o-l 

y(a  ) 

1.920 

1.357 

1.446 

Po(debyes) 

0.121 

0.166 

0.223 

In  order  to  estimate  |±q  and  y  for  systems  considered  here  it  is 
assumed  chat  y  ^  is  proportional  to  the  collision  diameter  ct  occurring 
in  the  interparticle  potential  and  that  is  proportional  to  both 
the  collision  diameter  and  the  number  of  electrons  in  the  outer  shells 
of  the  component  particles.  The  first  assumption  seems  reasonable, 
since  the  range  of  the  induced  moment  is  proportional  to  the  separation 
of  the  particles  when  their  electron  clouds  are  overlapping  and  hence 
to  the  collision  diameter  characteristic  of  the  interparticle  potential. 
The  second  assumption  also  seems  reasonable  since  the  magnitude  of  the 
dipole  moment  is  proportional  to  the  mean  separation  c?  charges  (which 
is  taken  to  be  y  and  the  mean  number  of  electrons  involved  in  the 
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collision  (taken  as  the  electrons  in  the  outer  shells  of  the  colliding 

tk 

particles).  The  inner  electrons  are  too  closely  bound  to  be  of  mucb 


consequence  to  the  induced  moment. 


Hence, 


and 


jonstant 


yc  ■  constant 


(33a) 


(33b) 


where  a  is  the  collision  diameter  and  N  is  the  number  of  electrons  in 
the  outer  shells  of  the  collidin  particles.  From  the  values  given  in 
the  Table,  for  example,  Eq.  (33a)  cairly  well  satisfied.  From 
Eqs.  (33a)  and  (33b) 


Oiv'1)B 


(34) 


where  the  subscripts  refer  to  two  different  mixtures.  For  the  three 
cases  given  in  the  Table: 


Since  the  induced  dipole  moment  is  taken  as  a  superposition  of 
two  terms  (the  rare-gas  atom  interacting  with  each  of  the  two  atomic 
components  of  the  molecule  respectively)  the  molecular  electron 
contribution  to  N  is  taken  as  one  half  of  the  number  of  outer  she1 l 
molecular  electrons  while  all  of  the  outer  shell  electrons  in  t 
rare-gas  atom  are  counted. 
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1)  A  -  He-Ne 


B  -  He-Ar 


JA 

—  -  0.73 

u 


0.71 


2)  A  ■  He-Ne 


B  -  Ne-Ar 


0.54 


cV1), 


0.47 


3)  A  -  He-Ar 


B  -  Ne-Ar 


0.74 


(ny-1)B 


0.68 


The  ratios  agree  fairly  well  for  these  cases. 

Values  of  y  and  for  H2*He,  Hj-Ar  and  Nj-Ar  systems  are  now 

estimated.  Equation  (33b)  Implies  y  /y  ■  a '/a..  Then 

A  B  BA 


^2 ’He  "  ^He-Ar^  aHe-Ar^CTH2-He 


1.6 


(35a) 


and 


VAt  -  <vHe.Ar)  V-Ar/0H2-Ar  '  1'J  *' 

VN0-Ar  "  ^VHe -Ar^  aHe-Ar^aM  -Ar  "  1,1  A 

*  l 


(35b) 


(35c) 
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to  two  significant  figures.  In  arriving  at  these  results  the  value  of 

the  collision  diameter  a  was  obtained  in  the  usual  way  as  the  arithmetic 

CIO!  0 

mean  of  the  values  for  the  component  gases.  '  (cr„  ■  2.46  A, 

**2  "He 

-Ar”  -Ar™  ^He-Ar"  ^.91  X.)  From  Eq.  (34)  , 


(Ny"1) 
-He  "  (Ny-1) 


H2-He 


He-Ar 


0.044  debyes 


He-Ar 


(36a) 


<W'\-Ar 

(.  )  SB  - —A -  (u J 

VAr  (My  )He.Ar  He-Ar 


0.16  debyes 


(36b) 


W'V-Ar 

(tO  - Zl -  <0 

N2-Ar  (Ny  )He.Ar  He-Ar 


0.27  debyes 


(36c) 


The  quantitative  results  for  the  absorption  spectrum  follow. 

ABSORPTION  SPECTRUM  FOR  Hg/RARE-GAS  ATOM 

First  consider  the  H2-He  system.  The  reduced  mass  m  for  H2_He  is 
approximately  4/3  amu  (atomic  mass  unit)  and  the  rotational  constant  B 
for  H^  is  59.3  cm  jn  this  example  and  the  following  the  tempera¬ 

ture  T  is  taken  to  be  of  the  order  of  room  temperature  (i.e.,  xT  =>*  200  cm 
Using  these  values  for  B  and  xT  it  is  easy  to  see  that  only  the  ground 
and  first  few  (about  three)  excited  rotational  states  of  H2  are  appreci¬ 
ably  excited  under  equilibrium  conditions.  The  nuclear  spin  of  H2  is 
hence  the  total  wave  function  must  be  antisymmotrical.  This  implies 
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that  the  relative  nuclear  statistical  weights  are  1  and  3  for  even 

and  odd  rotational  states  respectively.  The  polarizability  of  He  is 
-24  3 

—  0.21  X  10  cm  (Ref.  11)  and  the  quadrupole  moment  of  is 
2 

2s  0.49  ea  (Ref.  7)  where  e  is  the  electron  charge  and  a  is  the  Bohr 
o  o 

radius.  Substituting  numerical  values  in  Eq.  (32)  gives  (at 
kT  2:  200  cm*  *) 


aH  -He(k)  ~  1,53  X  10'46  [Q  +  I  S(J)]  n1n2<cm’1> 


(37) 


J=0 


where  Q  is  the  sum  given  in  Eq.  (32)  with  «*  J^,  -  0,1, 2, 3,  S(J) 

is  given  by  the  sunmand  in  Eq.  (32)  with  J  s  J^,  *  J^-2,  and 

the  wave  number  is  k  cm  *.  In  Fig.  1  the  contributions  of  the 
various  branches  and  the  sum  of  these  branches  are  plotted  as  a  function 
of  the  wave  number  k  (i.e.,  k  =  u)/2tt)  .  The  Q  branch  contains  contribu¬ 
tion  from  both  the  overlap  and  quadrupole- induced  dipcle  moments  (the 

overlap  term  being  ~  90  percent  of  the  Q  branch  contribution)  ,  while 

** 

only  the  quadrupole- induced  moment  contributes  to  S  branches. 

As  another  example  the  H0-Ar  system  is  considered.  The  reduced 

tm 

mass  H^-Ar  is  approximately  1,91  amu,  and  the  polarizability  of  Ar  is 
-24  3 

1.63  X  10  cm  (Pef.  12).  Substituting  numerical  values  into 
Eq.  (32)  gives  (at  kT  “■  200  cm  ^) 


In  all  cases  discussed  here  the  contribution  to  a  from  the  0 
branches  was  found  to  be  negligible  compared  to  the  contributions 
from  the  Q  and  S  branches. 

At  u>  =  E  T  (resonance)  Eq.  (32)  indicates  that  the  quadrupole 

Jr  i 

contribution  becomes  very  large.  I  —  log  (u)-E  ).  Strictly  speaking, 

q  Vi 


the  rotational  transition  lines  are  not  infinitely  short  but  have  very 
narrow  pressure  broadened  lines  («  1  cm"*)  at  room  temperature  and 
atmospheric  pressure.  The  absorption  at  resonance  is  finite;  I  ~  log 

*>  2  ^ 

[(tt/-ET  T  )  +  V  /4 3  where  y  is  the  half  width  of  the  line.  Due  to  the  lack 

Y  i 

of  knowledge  of  these  line  widths,  the  increase  in  the  absorption  spectrum 
very  near  a  resonance  is  indicated  schematically. 
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3 

c*H  _Ar00  —  3.75  X  10'46  [Q+  l  S(J)]  n1n2(ctn‘1)  (38) 

2  J=0 

In  Fig.  2  the  contribution  of  the  various  branches  and  che  sum  of  these 
branches  are  plotted  as  a  function  of  wave  number.  In  this  case  the 
absorption  coefficient  is  roughly  an  order  of  magnitude  larger  th<?.n 
the  H^-He  case.  This  is  due  mainly  to  the  larger  reduced  mass  of  the 
H2~Ar  system  and  the  larger  polarizability  and  induced  moment  of  Ar 
and  H2~Ar  system  as  compared  with  l^-He. 

ABSORPTION  SPECTRUM  FOR  t^-Ar 

Due  to  the  relatively  large  rotational  constant  of  (Bq  *  59.3  cm  *)  f 

the  absorption  spectrum,  at  room  temperature,  can  be  resolved  into 

two  parts;  a  pure  translational  component  (Q  branch),  and  a  rotation- 

translation  "omponent  (S  branch) .  On  the  other  hand,  in  tbi  case  of 

N  the  rotational  cor'tant  is  much  smaller  than  for  H  (B  2:  2  cm  *) , 

2  2  Kj 

with  the  resu.t  that  at  room  temperature  the  Q  and  S  branches  completely 

overlap.  That  is,  a  resolution  of  the  absorption  spectrum  of  ^/rare- 

gas  atom  systems  into  a  translational  component  and  a  rotation- 

translation  component  is  not  possible.  The  N0-Ar  system  is  considered. 

The  reduced  mass  of  Nj-Ar  is  approximately  10.4  amu,  thr  quadiupole 

2 

moment  of  N2  is  =“  1.1  ea^  (Ref,  13).  The  nuclear  spin  of  N2  is  1, 
which  implies  that  the  relative  nuclear  statistical  weights  are  2  and 
1  for  even  and  odd  rotational  states  r*  jpectively.  At  room  temperature 
it  is  necessary  to  include  the  ground  and  about  the  first  twenty 
excited  rotational  states  in  the  statistical  average.  Substituting 
numerical  values  into  Eq.  (32)  gives  (kT  _0O  cm"*) 


Arrows  along  abscissa  indicate  frequencies 
of  unperturbed  rotational  transitions 
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Absorption  coefficient  for  H2-Ar  system  at  room  tsmperature  with  contributions  of  Q  and  S  branches  indicated 
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<VM  a  (k)  a:  2.6  X  10  R(k)  n.n  (cm  1)  (39' 

Nj “Ar  i  2 

20 

where  R(k)  ■  Q  +  £  S(J).  The  absorption  coefficient  for  JL-Ar 

j-o 

is  plotted  in  Fig.  3  (the  S  branch  contributes  about  90  percent 
of  the  total  to  R(k)).  Very  close  (<<  1  cm  ^  to  a  transition 
line  the  absorption  becomes  large,  but  this  occurs  in  a  very 
narrow  region  about  each  line  (compared  to  the  spacing  between  them). 

Due  to  lack  of  knowledge  of  the  relevant  line  widths,  the  curve 
shown  in  Fig.  3  is  the  envelope  of  points  away  from  rotational  transi¬ 
tion  lines.  Very  close  to  the  resonances  the  absorption  spectrum 
has  a  series  of  very  narrow  spikes.  The  contribution  of  these  resonance 
lines  to  the  integrated  absorption  coefficient  is  negligible,  since 
these  lines  are  so  narrow;  therefore  they  are  not  indicated  in  Fig.  3. 

From  the  N^-Ar  spectrum  an  estimate  of  the  absorption  in  N^-N^ 
that  involves  single  rotational  transitions  may  be  easily  obtained. 
Single  (as  opposed  to  double)  rotational  transitions  are  those 
which  involve  a  rotational  transition  in  one  of  the  colliding  mole¬ 
cules  but  not  in  both.  To  obtain  this,  one  substitutes  the  polariz¬ 
ability,  reduced  mass,  the  collision  diameter  for  the  system, 

and  an  overall  factor  of  two  for  that  of  the  N^-Ar  system.  This 
method  treats  one  of  the  modules  essentially  as  a  rare-gas  atom — 
that  is,  as  having  no  internal  energy  levels  to  excite  it.  It  is 
found  that  at  the  peak  of  the  spectrum  (—  100  cm  ^ ,  room  temperature 
and  atmospheric  pressure)  the  absorption  coefficient  is  of  the 
order  of  0.8  X  10  ^  cm  \  Bosomworth  and  Gush^  report  an  experi¬ 
mental  absorption  coefficient  of  ~  4  X  10  ^  cm  ^  at  ~  100  cm 
Presumably  the  reason  for  this  discrepancy  of  values  lies, 


\ 


oo 


Fig.  3— Absorption  coefficient  for  N  -Ar  system  at  room  tempe 


33 


in  part,  in  the  failure  of  our  analysis  to  include  the  contribution 
of  double  rotational  transitions. 


ESTIMATES  FOR  N2-N2  ABS0RPTI0N  IN  THE  FAR  WINGS 

For  frequencies  oj  »  hT  and  E  the  momentum  transfer  Q  becomes 

JfJi 

large  (Qa  »  1  where  a  is  one-half  the  intemuclear  spacing  0.5  A 

for  N^) .  This  can  be  seen  as  follows:  Q  =  | k  k  j  and  k'  /2m  =  k  /2m 

2  u 

+  u>  -  E  ;  k  /2m  hT,  hence  for  uu  »  kT,  E  .  implies  Q  «*  (2muu)  . 

Vi  ,  Vi 

Thus  for  large  u,  Q  *.r,  At  the  CO^  laser  wavelength  (k  2;  10  cm  ) 

we  find  Qa  ~  1,  and  hence  the  overlap  moment  is  angle-dependent 

(i.e.,  e  1)  .  Large  Q  implies  close  collisions  and  hence  the 

quadrupole- induced  moments  do  not  contribute.  Also,  since  the  overlap 

moment  is  angle-dependent  the  ^*^2  system  no  longer  possesses  a  center 

of  symmetry  and  can  have  a  nonzero  overlap- induced  dipole  moment.  By 

a  process  similar  to  that  which  led  to  c,qs.  (35)  and  (36),  it  is 

o 

found  that  (o„  =  3.75  A) 

N2'N2 


and 


9-1 


VN2-N2  “  <VHe.Ar)  CTHe-Ar/aN2-N2  “  1,05  A 


(NY  }n2-n 

(u  )  »  - : -  (u  )  ■  0.11  debyes 

°  VN2  >He-Ar  °  He‘Ar 


(40) 


(41) 


Instead  of  using  Eq.  (12),  we  take  the  overlap-induced  iipole 
moment  of  the  two  molecule  system  to  be 


g!( r )  =  u^)  +  u^')  +  u(r2)  +  u(r2  ') 


(42) 
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The  vectors  ,  r^',  r^,  and  r^J  are  defined  in  Fig.  4.  Letting 
r  be  the  vector  between  the  centers  of  mass  of  the  two  molecules,  and 
calling  the  molecular  internuclear  separation  vector  2a^  and  28^, 
where  the  subscripts  1  and  2  refer  to  the  two  molecules  respectively, 
it  is  seen  that 


r  -  a„  -  r„  -  a, 
2  2  1 


r  +  a2  •  rl  ‘  *1 


?  -  a.  -  r„ #  +  a.  ■  0 

2  2  1 


— »  *4  — *  <  —4 

r  +  &2  -  -  0 


(43) 


The  integral  over  the  relative  coordinate  (r)  is  the  same  as  before, 
that  is,  in  place  of  Eq.  (17),  it  can  be  shown  that 


H  =/(JH)^1)(JM)^2)1uz1(JM)^1)(JM)^2)\  -  AqQ2 


f|e1*Sl|iN 


(1) 


x  <  f|eiQ‘a2|i)(2)  +  (f!eiQ'al|i)m 


'(1) 


(2) 


+  terms  a^  -a^,  a^  -»  -a^ 


(44) 


where  the  sub-  and  superscripts  1  and  2  refer  to  each  of  the  two  molecules 
respectively  and 

_3/2  9  9 

'ts _  e-(Q2/4‘,2) 


A  ■ 

o 


2 


(45) 
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nary  molecular  collisions 


I 


Taking  the  z-axis  parallel  to  the  polarization  vector  of  the  photon, 

iQ  •  g 

we  expand  e  in  terms  of  spherical  harmonics;  that  is. 


-  *n  f  iJ  £  ^(0.)  4‘V 


where  Y.  are  spherical  harmonics,  and  Cl  and  fL  are  the  angular 
JM  a  V 

coordinates  of  the  vectors  a  and  Q  respectively.  Because  of  the 

symmetry  of  the  molecule  about  its  center  of  mass,  only  even  J's 

■k 

will  contribute  (as  can  be  deduced  from  Eq.  (44)).  Also,  since 

Qa  ^  1  the  asymptotic  form  of  the  spherical  Bessel  function  is  used: 

j  (Qa)  ~  sin  (Qa  -  4^) /Qa  =  (-1)^“  sin  Qa/Qa  for  even  J. 

J  ^ 

The  average  over  the  azimuthal  quantum  numbers  takes  the  form 


(2Jj+l)  (1)  (2Ji+1)(2) 


(M.)  (M  )  (M  )  (M  ) 
11  2  1  2 


Substituting  Eq.  (46)  into  Eq.  (44)  and  obtaining  the  indicated 
rotational  matrix  elements  (which  can  be  expressed  in  terms  of  Clebsch- 
Gordan  coefficients) ,  the  operations  indicated  in  Eq.  (47)  are  then 
performed.  The  sums  over  the  M  values  can  be  obtained  by  using  the 
orthogonality  conditions  on  the  Clebsch-Gordan  coefficients,  and  it 


is  found  that 


2  2 
H  “  16Q 

z 


sin  Qa  2  sin  Qa  2  mj  (2) 

H?r]  f1i’° 


Since  J  is  even,  it  can  be  deduced  that  must  be  even, 
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where 


F(i,f) 

In  the  following, 
that 


Y  (2J+1)2  C2(J.,J,Jf,00) 
J(even) 

a.  Replacing  Eq.  (30),  it  is  found 


a(o))  ~  (2tt)  ^(4m/3ea^)  f  7^  sin^  Qa  e  ^  e  ^ 

J  0  ^ 


(49) 


"  (2tt)  ^(4m/6ea^)  eG^  ^2  f  J  sin^  Qa 

JO  Q2 


“p  ['q2(^  +  i>  ■  3] 


2y 


2 

-  (2TT)"4(4m/6ea4)  ecA  I^oj) 


where 


and 


uw '  r 

Jo 


dx 


>[x(A+B)  +  J] 

*  sin^Cx^ 


2  2  2 
A  =  A  /2y  =  itou/y 


B  =  /4  =  tu/4xT 

C  =  Aa  =  (2ma))^a 


For  uj  —  CO^  laser  frequency,  C  >  1  and  hence  the  sin^  term  in  the 
integrand  of  1^  oscillates  rapidly  over  the  region  where  the 


4Q 


(50) 


exponential  is  nonzero.  Hence,  its  average  value  may  be  used, 
yielding 


t 
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where 


3 

8 


CO 


dx 

X 


-[(A+B)x  +  J] 


■  4  <i+S>  V2"  ^ 


5  -  B/A  «  4mxT/v2 


(51) 


Hence  it  is  found  for  03  »  xT,  E  .  that 

JfJi 

,  2 2  oj/2kT  „  |  os  /  4mxT  | 

(«/2)«onoe  h\paJl+-T) 


*VN2 


(a))  3: 


8(2n»KT)Va4  (l  +  ^f) 


20 


-  ?=  J4(J4+1) 


Z_1  ^  6  HT  ^  1  F(i.f)gi(2Ji+l)[  (52) 


J^O 

Jf-Ji+Aj 


where 


20  -  J.(J.+1) 

Y  e  hT  gi(2Ji+l) 


V0 


-3  * 

F(i,f)  is  given  by  Eq.  (49),  and  the  molecular  density  is  n  cm  .  In 
obtaining  Eq.  (52)  the  term  (2J^+l)e  has  been  dropped  since 

03  »  xT.  For  systems  the  reduced  mass  is  7  amu  and  a  *  .505  A 

Hence  from  Eqs.  (40),  (41)  and  (52)  it  is  found  that  at  room 
temperature 


a  o  2 

The  factor  n  /2  appears  in  Eq.  (52)  (instead  of  n  )  in  order  to 

count  each  collision  once. 
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*N  -N  (k)  “  4,5  X  l°'U2  ^o(3-17k)ekMj.  (cm-1)  (53) 

~  3.2  X  10  "42  n2  k5  e"2,92k,  for  k  »  1/3.17 

—  .o 

where  k  ■  10  k.  In  obtaining  these  results  only  those  final  states 
corresponding  to  AJ  ■  0,  +  2  have  been  included.  Since  the  overlap- 
induced  moment  falls  off  very  rapidly  with  increasing  interparticle 
separation  (the  moment  having  a  maximum  for  r  =  y  1  molecular 
separation),  angular  momentum  considerations  imply  that  AJ  =  ~  2 

for  .  This  can  be  seen  by  estimating  the  order  of  magnitude  of  the 
relative  angular  momentum  of  the  colliding  particles,  mvr,  with 
v  ~  (nT/m)^  r  ~  y  and  noting  that  the  photon  contributes  one  unit 
of  angular  momentum.  At  the  CO^  laser  wave  number  k  Ss  9.44  it  is 
found  that 

°^kC0  )  -3'3X  10"54n2  (cm'1)  (54) 

-  2.5  X  10‘15  cm'1 

19 

at  atmospheric  pressure  (n  ~  2.7  X  10  ),  which  is  very  small.  Even 

at  high  pressures  the  absorption  coefficient  remains  small.  For  example, 
at  100  atmospheres  the  absorption  coefficient  at  10. 6^  is  —  2  X  10  11  cm  1, 
Finally  it  is  noted  that  the  absorption  coefficient  is  relatively 
insensitive  to  temperature;  therefore  no  large  change  in  a  is  obtained 
by  varying  the  temperature. 
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V.  CONCLUSIONS 

The  results  obtained  for  Nj/rare-gas  atom  and  absorption 

in  the  far  wing  were  deduced  from  a  model  which  assumes  a  specific 

functional  form  for  the  overlap-induced  dipole  moment.  The  dependence 

(4) 

on  interparticle  separation  is  of  the  form  used  by  Levine,  while 
the  angular  dependence  is  introduced  by  assuming  that  the  total  over¬ 
lap  moment  is  a  superposition  of  one  particle  interacting  with  each 
of  the  atomic  components  of  the  molecule  respectively,  the  latter 
assumption  having  important  effects  for  close  collision,  when  details 
of  the  intermolecular  potential  become  important.  For  close  collisions 
the  assumption  of  plane  waves  to  describe  the  relative  translational 

motion  of  the  two  particles  is  offset  by  assuming  that  the  dipole 

.  2  2 

moment  goes  to  sero  for  small  r,  i.e.,  p,  =  M>oYre  ^  r  •  The  factor  yr 

is  used  to  simulate  the  effect  of  the  scattering  wave  functions. 

(4) 

Levine  has  shown  that  this  model  is  relatively  insensitive  to  the 
effects  of  the  interp.. rticle  potential  for  wave  numbers  near  the  peaks 
in  the  rare-gas  spectrum  (100-400  cm  *) .  No  experimental  results  as 
yet  exist  for  the  far  wings,  and  hence  the  accuracy  of  an  extrapolation 
of  the  overlap  model  into  the  far  wings  cannot  be  ascertained.  The 
results  obtained  here  for  Nj-Nj  absorption  at  k  °*  1000  cm  *  have  Leer 
based  on  an  extrapolation  of  this  model.  The  absorption  coefficient 
obtained  is  very  small  ar^  it  is  not  f'lt  that  any  refinement  of  the 
model  will  produce  orders  of  magnitude  changes  in  a.  Therefore  it  is 
concluded  that  this  mechanism  does  not  play  an  important  role  in  the 
atmosphere  for  the  absorption  of  10. 6n  radiation  by  nitrogen. 
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Appendix 


EVALUATION  OF  I(cu) 


Consider  the  integral  I(»)  defined  in  Eq.  (28): 


-Ax  -  * 

,  **  x  2k 

dx  e  sin  Cx 


A,  B  and  C  are  given  in  Eq.  (28)  At  x  -  0  and  x  -•  «  t.ie  integrand 

approaches  zero  exponentially  and  has  a  maximum  for  x  =  (B/A)^.  Upon 

examining  the  numerical  values  of  parame.ers  for  the  cases  discussed 

2 

in  the  text  it  can  be  seen  that  the  sin  Cx  term  is  a  slowly  varying 
function  of  x  it.  the  region  where  the  exponential  term  is  nonzero. 
Hence,  upon  substituting  those  values  of  x  in  which  the  exponential 


term  is  maximum  into  the  sin  term. 


I(o>)  ="  sin2[C(B/A)*]  f  dx  e 

J  n 


‘Ax  -  l 


'  -  »l":CC(l/A)t]  jj_  {Jq  *  e  *}  (56) 


The  formula 


2  sin2[C(B/A)*]  |j  {k0[2(AB)^]]- 


.00  ,  -Ax  - 


f  * 

Jft  x 


-  2K0C2(AB)%] 


has  been  uced.  *  It  can  be  shown  that 


aK0[2(AB)^]  ,  , 

- — - -  -  (B/ A)*  K1[2(AB) H 


Hence,  from  Eqs.  (56)  and  (57), 


Ktu)  a:  2  sin2[C(B/A)^]  (B/A)*  ^^(AB)*] 


(58) 
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